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Abstract We report on the design of a segmented lin- 
ear Paul trap for optical clock applications using trapped 
ion Coulomb crystals. For an optical clock with an im- 
proved short-term stability and a fractional frequency 
uncertainty of 10^^^, we propose ^^^In+ ions sympa- 
thetically cooled by ^^^Yb^. We discuss the systematic 
frequency shifts of such a frequency standard. In partic- 
ular, we elaborate on high precision calculations of the 
electric radiofrequency field of the ion trap using the fi- 
nite element method. These calculations are used to find 
a scalable design with minimized excess micromotion of 
the ions at a level at which the corresponding second- 
order Doppler shift contributes less than 10^^® to the 
relative uncertainty of the frequency standard. 



1 Introduction 

Since the first demonstration of a single trapped ion 
in 1980, experiments with trapped single particles have 
led to some of the finest spectroscopic measurements in 
physics [UllllS] and radiofrequency (rf) Paul traps have 
become outstanding working tools with unsurpassed ac- 
curacy in frequency metrology [HIS]. 

The location of the single ion and thus, electro-magnetic 
fields sensed by it, can be controlled and measured at the 
nm scale, leading to a superior precision of optical ion 
clocks with a fractional systematic frequency uncertainty 
now reaching down into the 10~^* regime Still, the 
low signal-to-noise ratio of the single ion interrogation 
limits the obtainable short-term stability and puts up 
high demands for the stability of the optical clock laser. 
One approach to improve the short-term stability of ion 
clocks is to investigate narrower atomic transitions with 
mHz linewidth and thus, higher quality factor. Here, for 
Fourier limited spectroscopy a short-term linewidth of 
the clock laser at the mHz level is required. Due to the 
limited quantum information obtained in single ion spec- 
troscopy, long integration times are necessary to lock the 



laser onto the atomic signal. For example, resolving and 
locking a laser to a mHz wide clock transition by stan- 
dard quantum jump spectroscopy, would require a clock 
laser frequency stability in the low 10^^'' range over sev- 
eral minutes, in order to reach the quantum projection 
noise limit of the single ion [6 . A major effort is made in 
the metrology community to push the stability limits of 
optical reference systems [HIHIIS] and faster, phase sensi- 
tive detections schemes could be thought of. Still, today's 
state-of-the-art clock laser stability and detection limit 
the achieved short term stability such, that integration 
times of tens of days to weeks would be necessary to 
reach a frequency resolution of 1 x 10~^^ with an optical 
ion clock. This high frequency resolution is required for 
improved fundamental tests of physics and applications 
in geodesy [TO l Ill ) . As for white laser noise, the stability 
of the frequency measurement averages down with the 
square root of both atom number N and integration time 
T pjl , increasing the number of ions for example to only 
10 could shorten the integration time already by a full 
order of magnitude. Together with improved laser stabil- 
ity and possibly new detection methods, increasing the 
number of ions within a well controlled ensemble of ions 
would substantially improve optical frequency standards 
and open up ways for new applications in navigation and 
geodesy [lO] . 

A pivotal question is, whether a larger sample of ions 
can be controlled sufficiently well to reach fractional fre- 
quency inaccuracies or long-term stabilities as low as 
10~^^ or below. So far frequency standards based on mi- 
crowave transitions in a cloud of buffer gas cooled ions 
have been investigated for a frequency standard with a 
long-term stability of 10~^^ [13] • A ring of ions in an 
anharmonic linear Paul trap has been proposed [14] for 
an optical frequency standard with improved short-term 
stability, conserving the long-term stability at the level 
of 10-15. 

For an optical ion clock at the accuracy level of 10^^^, 
the choice of atomic elements is limited as many ions 
have an electronic quadrupole moment 9 in the excited 
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state of the clock transition, which makes them sensi- 
tive to electric field gradients. In a linear ion trap with 
its large static electric field gradient and the Coulomb 
fields of neighboring ions, quadrupole shifts are on the 
order of 10 Hz. An exception are two-electron systems, 
where transitions between ^Sq and ^Pq states with 6 — 
are available. The most advanced ion clock today 5j 
uses such a system, ^''Al+, with 9 = 0. Here, the clock 
ion is stored in a linear Paul trap together with a sec- 
ond species with easily accessible atomic transitions. The 
clock signal is read out via the second ion in a quantum 
logic step ^15iil6] . In general, the field of ion trapping has 
been pushed tremendously by the developments in quan- 
tum computation with cold ions in the past few years. 
A high level control of up to 14 ions has been demon- 
strated in entangled strings of ions [T71ITB] , cooling to the 
three-dimensional ground state was implemented [19] 
and new scalable ion trap structures are being developed 
for larger quantum registers [2U[|21) . 

The most striking remaining problems for optical 
clocks are heating of the ion due to its interaction with 
fluctuating patch potentials |22j and driven motion (mi- 
cromotion '23') due to residual rf fields of the trap, both 
giving rise to significant second-order Doppler shifts. In 
the NIST Al/Mg ion trap a fractional frequency shift 
due to micromotion of 3 x 10~^^ was observed, when the 
ion was shifted along the axial trap direction by only 
3 /Ltm [S]. This excess micromotion (EMM) arises be- 
cause of finite size effects, misalignment and machining 
tolerances in the ion trap assembly, see section 4, and 
poses severe problems, when scaling up an ion optical 
frequency standard to many ions. 

Inspired by scalable designs of ion traps for quantum 
computation [50] we present a design study for an optical 
frequency standard based on linear chains of ions with 
a potential fractional frequency uncertainty of 10~^^ or 
below. In our scalable, segmented four-layer Paul trap, 
chains of up to 10 ions can be trapped in the Lamb- 
Dicke regime [24] in each trapping segment. It provides 
both, all degrees of freedom to control micromotion as 
well as optical access in three dimensions. In such a trap 
we plan to trap multiple chains of ^^^In+ ions. The in- 
dium ion clock transition can be detected directly us- 
ing quantum jump spectroscopy. In order to ensure ef- 
ficient trapping and cooling of the ionic crystals, addi- 
tional sympathetic cooling with ^^^Yb+ will be provided. 
The high precision finite element method (FEM) calcu- 
lations presented in this paper are used to develop an 
ion trap made of AIN ceramics wafers in our lab. In a 
prototype of this trap, based on a ceramic filled, glass 
reinforced thermoset (Rogers4530B'^'^), we have already 
successfully trapped chains of ^^^Yb+ ions for first tests 
on this design [2F. Clearly, there are many challenges 
ahead in building an optical clock based on Coulomb 
crystals with systematic frequency shifts controlled at 
the level of 10^^^ and issues of controlling the ion dy- 
namics, collisions with background gas, rf phase shifts, 



magnetic fields and optical imaging will be addressed 
in our experiment with mixed chains of ytterbium and 
indium ions in our scalable ion trap. 

In summary, in this paper we propose a novel opti- 
cal frequency standard based on linear crystals of ^^^In"*" 
ions sympathetically cooled by ^^^Yb"*". We focus on the 
design of the ion trap for such an optical clock with many 
ions. In particular, we address the problem of residual 
micromotion and required machining tolerances. As mi- 
cromotion causes extra heating in laser-cooled systems, 
this study is of general interest for high precision exper- 
iments, for example for investigating ion-atom interac- 
tions in degenerate quantum gases |26| . 

In section 2, we detail our approach to an optical in- 
dium frequency standard with ionic crystals and show 
that a long-term frequency stability of 10~^* or below 
with systematic frequency shifts at the mHz level is pos- 
sible. Section 3 gives general requirements and consider- 
ations regarding the ion trap design. Section 4 describes 
the high precision FEM calculations of the electric field 
of a segmented linear ion trap. The influence of the ge- 
ometry is discussed and critical tolerances are detailed. 



2 Systematic shifts in an optical indium 
multi-ion clock 

The single ^^^In+ ion is a well known, previously investi- 
gated candidate for an ultra-stable and accurate optical 
clock 27,28,29 . Its narrow clock transition ^Sq ^^Pq 
at 236.5 nm, with a natural linewidth of 7 = 0.8 Hz and 
electronic quadrupole moment 9 = 0, see flgurejl] makes 
it an interesting candidate for a scalable optical clock 
with many ions. The possibility to detect the quantum 
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Figure 1: Simplifled level scheme of ^^'^In+ and ^^^Yb+, 
showing relevant cooling, repump and clock transitions. 



information of the clock excitation directly via the '^i-'i 
state with a natural lifetime of 0.44 /is (7 ~ 360 kHz) [37] 
can facilitate the atomic signal read-out of a larger chain 
of ions, without the need of quantum logic technique [15) . 
Its transition wavelength of 230.5 nm can be generated 
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Table 1: Systematic relative frequency shifts in ^^^In+. 
For the uncertainty of the BBR shift the theoretical un- 
certainty is given. (*) a heating rate of 50 phonons per 
second at Vscc = 1 MHz is assumed, see text. 



Type of shift 


Shift \A,y\ 
in 10"^** 


Uncertainty 
in 10"^* 


Thermal motion (T=0.5mK) 


1 


1 


Heating* 


< 0.1 


< 0.1 


Excess micromotion 


« 1 


< 1 


2nd order Zeeman 


< 0.1 


< 0.1 


BBR at 300K 


13.6 


1 


DC Stark 


0.08 


0.08 



with standard diode laser technology and second har- 
monic generation [30j . Due to the lower cooling power 
obtained by laser cooling on this narrow transition and 
the increased rf heating in Coulomb crystals, it is still 
advisable to sympathetically cool the indium ions with 
a second ion species with a stronger cooling transition. 
In our case we choose ^''^Yb+, because of its long life- 
time in ion traps, the all-diode-laser based easy accessi- 
ble wavelengths and the similar ion mass, that allows for 
efficient sympathetic cooling. YbH"*" formed in collisions 
with background gas can be dissociated with the Doppler 
coohng light jSlj- In addition, the ^^^Yb+ isotope has no 
hyperfine structure and cooling is facilitated. Besides the 
enhanced cooling efficiency and control of ions, the pres- 
ence of a second species allows for the possibility to sym- 
pathetically cool during the clock interrogation in case 
of excessive heating, additional characterization of the 
trap environment, such as magnetic fields, and an alter- 
native clock read-out via quantum logic for comparative 
studies. 

Besides the absence of electronic quadrupole moment, 
the ^5*0 ~^'^Pa transition in ^^^In+ has the advantage of 
a very low sensitivity to environmental effects, which are 
summarized in table[T] In particular, indium profits from 
its heavy mass when considering relativistic frequency 
shifts (second-order Doppler shift) due to time dilation 
Avfd/i^ = — £'kin/'7ic^, where -Ekin is the kinetic energy 
of the ion, m its rest mass, c the speed of light and the 
frequency of the atomic transition. At the Doppler cool- 
ing limit Td = 0.5 mK of ^^^Yb+ this frequency shift 
amounts to Ai'td/t^ = —Ekin/'mc^ = — 1 x 10^^^, where 
Ekin = (5/2)fcBr is the kinetic energy due to ther- 
mal secular and micromotion in a linear ion trap |23| . 
If relative uncertainties below 10~^® are targeted, fur- 
ther cooling such as sideband cooling of the Yb"*" ion 
or direct cooling on the intercombination line of the 
ii5jj-^+ jQjj^ |34j can be implemented, where temperatures 
of 100 /xK can be reached. An additional kinetic energy 
can be brought into the system by excess heating due 
to thermally activated patch potentials and contamina- 
tions on the electrodes Assuming an enhanced 
electronic field noise with a power spectral density of 



Se = 10^^^ (V/m)^IIz~^ in our ion trap with electrode- 
to-ion distance of 0.7 mm, see figure 6 in |36) . we obtain 
a maximum relative frequency shift of 1 x 10""'^^ due to 
heating of the secular motion during the clock interroga- 
tion, corresponding to 5 phonons at an eigenfrequency 
of the secular motion 1^300= 1 MHz. Here we assumed 
a maximum clock interrogation time of 100 ms, limited 
by the natural lifetime of the ^Pq state of ^^^In+. Com- 
monly observed heating rates in such ion traps can be 
a factor 10 to 100 lower [3S|, giving relative frequency 
shifts of 10~^° and below. 

Regarding ac Stark shifts due to blackbody radia- 
tion (BBR), ^^^In+ has one of the smallest sensitivities 
among optical clock candidates. With a differential static 
polarizabihty of Aao — 5x 10^*Hz/(V/m)^ between the 
^So and ^Po state [33, it is comparable to ^'^Al+ [321133]. 
Here, dynamic corrections due to the frequency distribu- 
tion are neglected, as they are typically of the order of 
a few percent or below. At room temperature this cor- 
responds to a blackbody shift of the clock transition of 
1.36(10) X 10~^^. For example, determining the temper- 
ature of the trap environment to T = 30 ± 10 °C will 
give a relative clock uncertainty due to the blackbody 
shift of 2 X 10~i^ 

Two major systematic frequency shifts arising from 
imperfections in the rf trapping potential are excess mi- 
cromotion (EMM) and dc Stark shifts. If the driving 
radiofrequency field E'j.f of the trap cannot be zeroed 
at the place of the ion, it amounts to a second-order 
Doppler shift of Avtd/i^ = —e^E^^/ic^m'^n^f, where e is 
the electronic charge, E^f the amplitude of the electric 
rf field at the place of the ion and i7if /27r the trap drive 
frequency 23J. This effect is the dominating uncertainty 
in today's best ion clocks and aggravates for a larger 
number of ions. Section 4 of this paper addresses the 
influence of finite size effects of the trap and imperfec- 
tions in the electrode machining on the rf field in detail. 
The targeted upper limit in our trap design for resid- 
ual rf fields on the trap axis is a maximum of 90 V/m, 
which corresponds to Ai/td/i' = — 1 x 10^^^ for a trapped 
ii5jj^-h jQj^ ^].g^p frequency i7i.f/27r = 20 MHz, as dis- 
cussed in section 3.1. In addition, any residual electric 
field seen by the ion, will give rise to a dc Stark shift 
I'g in the clock transition. Here, ^^^In+ profits from its 
low static differential polarizabihty Aao t32\ At an rf 
field amplitude of 90 V/m this Stark shift corresponds 
to vs = ^ao (E^f) = 8 X 10-20 x u. 

Lastly, we consider the influence of static and dy- 
namic magnetic fields in the ion trap. The linear Zee- 
man shift from the = 9/2 to tof = 7/2 states 
amounts to 6.36 kHz/mT |38j and is due to hyperfine 
mixing of ^Pq and ^Pi states (the nuclear spin of ^^^In 
is / = 9/2). It can be measured and subtracted by al- 
ternatively pumping the ion into the stretched states of 
the ground state with opposite magnetic moments |37l 
13]. The second-order Zeeman shift is given by Av = 
f3(^B'^^, where for the alkaline-earth like system /3 = 
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2/i^/3/i^Z\fs — 4-1 Hz/mT^, where [1b is the Bohr mag- 
neton and h the Planck constant. Owing to the large 
fine-structure splitting Z\fs — 3.2 x 10^"^ Hz of the ex- 
cited triplet states, indium has an advantageously low 
second-order B-field dependency, 10000 times lower than 
171YI3+ and 20 times lower than ^^Al+. Due to unbal- 
anced currents in the ion trap alternating magnetic fields 
with B^^g = 2.2 X lO-^^T^ have been observed m ion 
traps [51. For ^^Al+ this accounts for an ac Zeeman 
shift of Avjv = 1.4 X 10"^^. For the "5ln+ ion with 
Avjv = 7 X lO"^'^ this effect is negligible at the level 
of 10~^*. While static and dynamic B-fields can be de- 
termined quite accurately and frequency shifts can be 
taken account of, for an optical clock based on many 
ions, requirements on the homogeneity of magnetic fields 
become important. With a linear Zeeman shift of 6.36 
kHz/mT variations in magnetic field amplitude should 
be less than 16 nT across the ion chain to ensure that 
the broadening of the atomic line is less than 0.1 Hz. 
Distortions of the line profile due to spatially varying 
systematic shifts will have to be evaluated carefully to 
avoid locking offsets. The fast interrogation cycle, that 
is possible with a multi-ion frequency standard, will be 
advantageous when both stretched states are probed, in 
order to avoid offsets due to temporally varying mag- 
netic fields. 

The above considerations show that it should be pos- 
sible to evaluate the clock frequency of a larger sample 
of ^^^In+ ions with a fractional frequency uncertainty 
of 10~^*, assuming that a sufficiently ideal trap can be 
machined, which will be addressed in more detail in the 
following sections of this paper. We have not addressed 
other potential systematic frequency shifts like Doppler 
effects of first-order [3] or phase-chirps in acousto-optic 
modulators, which are of general technical origin and 
need to be evaluated in each individual experimental 
setup. Frequency shifts due to collisions with the back- 
ground gas have to be studied, but are typically in the 
10~^^ range when operating the ion clock at a back- 
ground pressure in the low 10^^ Pa range [Jj. 

3 Trap design 

For an optical clock with low systematic frequency shifts 
and competitive long term stability it is important to 
have a high level control of the electrostatic and mag- 
netic fields seen by the ions over the whole trapping 
region. This includes both necessary external control 
and 3D optical access to measure micromotion and field 
shifts. We consider a scalable segmented trap design, see 
fig. [2] in which ions can be trapped either in a long linear 
string with an axial anharmonic dc potential [39^ or, in 
an harmonic axial confinement, in separate "buckets" of 
up to 10 ions in every other trap segment, that can then 
be controlled independently. Compared to linear ion trap 
designs with separate dc electrodes on the trap axis for 



the axial confinement [40 , segmented rf electrodes [4TJ 
W2\ give a more ideal behavior of the rf field lines along 
the axial trap direction and allow for a high precision in 
the alignment of trap electrodes, if machined from one 
part. Currently a laser cutting and metallization process 
for low rf loss ceramic wafers is being developed in our 
lab. Four of these wafers can be precision aligned and 
stacked on top of each other to give the trap assembly. 

Figure 1(b) shows the top view onto one pair of rf 
quadrupole electrodes, that can be laser cut from one 
thin wafer. To the segmented rf ground electrodes differ- 
ent dc potentials can be applied for the axial confinement 
as well as for static stray field compensation in the ax- 
ial and one radial direction. To obtain both degrees of 
freedom to compensate for residual micromotion in the 
radial plane, our design includes an additional electrode 
wafer with compensation electrodes that overlie each rf 
ground/dc electrodes with a small spacing in the range 
of 0.1 mm to 0.25 mm, see fig. 1(a). The dc field gener- 
ated by these extra compensation electrodes as well as 
their effect on the rf field is described in section IITTI Our 
trap design is rather open and offers optical access in 
the (a;, z) plane and in the (y, z) plane for laser beams 
necessary for laser cooling, exciting repumping transi- 
tions, state manipulation and probing. In x direction, 
with a typical electrode thickness between 0.2 mm 
and 0.4 mm, a solid angle fraction larger than 4 % is 
available for optical fluorescence detection. 

3.1 Trap parameters 

In order to reduce heating of the secular motion of the 
ions we have chosen a rather large distance d = 0.7 mm 
from the electrodes to the trap center, since observed 
excess heating rates scale as d~'^ The geometry 

of the linear quadrupole electrode assembly has unity 
aspect ratio, with a distance of 1 mm between the tips 
of opposite electrodes, as shown in flgure[2j[a). For this 
case, the quadrupole term of the radial rf potential is 
only a factor 1.3 lower than in an ideal Paul trap with 
hyperbolically shaped rf electrodes. In this sense, the 
trap geometry is relatively efficient and sufficiently large 
trap frequencies of secular motion in radial direction Vr 
can be realized despite the large distance to the trap 
electrodes. With an rf voltage amplitude of 1.5 kV at 
an rf drive frequency Q-d/ii: — 20 MHz we obtain trap 
frequencies for the secular motion of v^.Yh = 1-25 MHz 
and Vr.in = 1-75 MHz for ^'''^Yh+ and "5ln+ ions, re- 
spectively. 

An important criterion for first-order Doppler free 
clock spectroscopy is the Lamb-Dicke condition with 7/ = 
kAx < 1 [24]. Here, Ax is the rms spread of the position 
of the ion and k is the modulus of the wave vector of the 
light exciting the transition. At an ion temperature cor- 
responding to the Doppler limit for Yb+, Td =0.5 mK, 
the Lamb-Dicke parameter in radial direction takes a 
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Figure 2: Schematic drawing of our trap design based on 
four segmented electrode wafers: (a) side view; (b) top 
view onto one inner wafer with a rf trap electrode pair. 
The following lengths of the trap geometry are defined: 
electrode thickness tc\ total length of the trap electrodes 
it; segment length Is, and width of isolation slit in the 
rf ground / dc electrodes w-^ which is identical to the 
width of the corresponding notches in the rf electrodes. 



value of ?7 = kAx = 0.45 for the indium clock transition 
at 236.5 nm. For the ground state of the harmonical os- 
cillator, with Vr^in ~ 1-75 MHz, 77 — 0.13. 

The confinement in axial direction is generated by 
applying a dc potential between neighboring segments in 
the rf ground electrodes. For an harmonic confinement 
inside the Lamb-Dicke regime, every other segment can 
be used as an ion trap with one segment for dc con- 
finement in between. We find that the curvature of the 
axial trapping potentials has a shallow maximum as a 
function of segment length Ig close to 1 mm dropping 
off by 10% at 0.8 mm and by 7% at 1.2 mm. Regard- 
ing the number of traps obtained per unit length of the 
quadrupole electrodes it seems desirable to use a smaller 
segment length. However, for segment lengths Is smaller 
than 1 mm neighboring traps become less well isolated 
from each other in the sense that varying compensation 
voltages in one trap are sensed by neighboring ion traps. 
In addition, the depth of the static potential well on the 
trap axis decreases monotonically with decreasing seg- 



ment length Is- As will be described in section 4.3 the 



choice of the segment length also influences the axial 
component of the rf field, as far as the contribution from 
the notches and isolation slits in the quadrupole elec- 
trodes is concerned. For our design, we consider segment 
lengths in the range from 1 mm to 2 mm. 

In an ideal linear Paul trap the rf field vanishes along 
the symmetry axis of the ion trap and driven micromo- 
tion of the ion can be minimized for a linear string of 
ions. To prevent transitions of the linear chain of ions to 
zig-zag, helical and more complex crystal structures, the 
ratio of radial to axial trap frequencies of secular mo- 
tion Vr/vz has to be kept above a critical value, which 
depends on the number of ions N in the crystal. The var- 
ious phase transitions in the ion Coulomb crystal were 
investigated numerically and analytically [47 , 45 , 46 . Es- 
timates for the critical ratio v^/vz for the transition from 
the linear crystal as a function of N from different au- 
thors are in agreement with each other and are compiled 
in reference [35]. Using the criterion Vr/vz > 0.73A^*''^^ 
obtained by Steane |44| . we find, for example, Vj-jvz > 3 
for five ions and Vj-jvz > 5.3 for ten ions. Choosing 
1.75 MHz for Vr the linear confinement of ten ions is 
possible for Vz < 0.33 MHz. 



3.2 Design criterion: micromotion in linear rf trap 

A principal limitation in today's optical ion clocks are 
frequency shifts induced by excess micromotion |23| . The 
observed EMM can originate from uncompensated static 
stray fields that shift the ion off the nodal line of the rf 
field, residual phase shifts of the rf potential on oppos- 
ing electrodes and imperfections in the ion trap geome- 
try that lead to a non-vanishing rf field along the trap 
axis. Our design allows to compensate static stray fields 
for all degrees of freedom in the radial plane to zero 
micromotion radially. RF phase shifts can be suppressed 
by carefully designed electronic circuits, symmetric leads 
and low resistivities. The electronic circuit scheme of our 
scalable trap design is currently tested on a first proto- 
type trap \25\ . In this paper we focus on the minimiza- 
tion of residual rf fields i?if.z along the trap axis, that 
are due to the trap geometry and cannot be zeroed, thus 
severely limiting the available trapping region, where mi- 
cromotion is negligible. 

The dominating frequency shift due to EMM is the 
second-order Doppler shift, while dc Stark shifts induced 
by residual rf fields are roughly an order of magnitude 
smaller in the case of ions with low static polarizabilites, 
such as ^^^In+ and ^''Al+, see section 2. We therefore 
assume as a trap design criterion a maximum second- 
order Doppler shift of lAiyfd/i^l = 10~^* due to residual 
axial rf fields. 

As already pointed out in section 2, this relativistic 
frequency shift contributes with Avtd/'^ = ^^'kin/l^^c^) 
to the fractional uncertainty of the optical clock. Here, 
the kinetic energy Skin = ^mT^ + -^kin™^ 
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mass m is given by the contributions E^^^ from the 
secular motion and E^^"^'^ from the micromotion of the 
ion. Assuming the secular motion to be thermalized at a 
temperature T one obtains E^^^'^^ = |A:bT for the three 
degrees of freedom of motion with ks denoting the Boltz- 
mann constant. For the radial degrees of freedom in a 
linear Paul trap, assuming the confinement is given by 
a purely quadratic pseudopotential 

tf' = e2^y(4mf?2^), (1) 

the thermal contribution of micromotion to the kinetic 
energy is equal to the kinetic energy contribution from 
the secular motion i?^™™^ = ffceT, a result obtained 
from both classical and quantum-mechanical treatment 
of the ion motion [48] . This is the case when the EMM 
is reduced to a negligibly small level in the radial direc- 
tions and when the contribution of the static potential to 
the radial confinement can be neglected. The static po- 
tential has a weakening effect on the radial confinement 
and thus results in an increase of the micromotion con- 
tribution to the kinetic energy which can be calculated 
using equation 10 in reference [53]. In an ideal trap the 
axial confinement is purely given by the static potential. 
The micromotion contribution to the kinetic energy in 
axial direction can only stem from excess micromotion 
due to a residual axial component of the rf field i^rf.z 
along the trap axis. Neglecting the contribution of the 
static potential to the radial confinement we obtain 

u 2toc2 4m2 02/22' 

when the excess micromotion in radial directions is re- 
duced to a negligibly small level, which is possible using 
techniques described by Berkeland et al. [IS] and taking 
care in the design of the rf circuit to minimize rf phase 
shifts between the quadrupole electrodes to a negligibly 
small level. 

As the second term on the righthand side of equa- 
tion [2] suggests, the axial component of the rf field -Erf.z 
has to be considered separately. Although a single ion or 
a linear chain of ions may also be displaced in axial direc- 
tion in the trap by applying a dc electric field in search 
for a minimum in |i?rf,z|i the existence of a minimum and 
whether the minimal magnitude of -Erf^z is small enough 
is determined almost entirely by the design of the trap 
and the strict keeping of tolerances in its construction. 
Hence, in the design of the trap we look for geometries 
which minimize the axial component of the rf field and 
analyze the sensitivity of different sizes in the geometry 
in order to state the tolerances for the construction. 

For indium ions using an rf drive frequency fl^ijli: = 
20 MHz and an rf amplitude of 1.5 kV a residual axial rf 
field amplitude |£'if,z| of 90 V/m amounts to a relative 
frequency shift of Avtdjv = — 1 x 10^^^ due to non- 
thermal EMM, second term in eq.[2] In order to minimize 
the axial micromotion for a linear chain of 5 to 10 ions 



to that level, in each trapping segment a region of a 
length of 50 /im to 100 /im is required along the z axis, 
where the magnitude of the axial rf field is no larger than 
90 V/m. 

4 Trap field calculations 

An ideal infinitely long linear Paul trap has no axial 
rf field component. In segmented linear Paul traps of 
finite length used in practice, a nonzero axial rf field 
component |i?rf,z| on the trap axis appears for several 
practical or technical reasons, which we investigated in 
order to minimize their effect. These are: the finite length 
of the quadrupole electrodes; insulation gaps in the dc 
electrodes; deviation from the ideally parallel alignment 
of the quadrupole electrodes to the trap axis; differences 
in the geometry of the electrodes and insulation gaps 
that breach the symmetry. 

For the trap field calculations in sections [4.2| through 
|4.6| an electrode geometry similar to figure [2] is used, 
with a varying number of segments, as will be indicated 
in the text. Only the layers of extra compensation elec- 
trodes are omitted, in order to preserve symmetries and 
to separate the different effects. The effect of these extra 
layers on the rf field is discussed in section [477] 

In order to obtain the ponderomotive potential given 
in equation 1, the electric field E^i of the linear Paul 
trap is calculated. The Laplace equation is solved for 
boundary conditions given by the electrodes. As this is 
a linear homogeneous differential equation, a solution 
for a specific boundary condition can be written as a 
linear combination of solutions of this geometry. It is 
worthwhile noting that a better understanding of the 
effects can be gained by decomposing the electrostatic 
potential problem linearly into one component contain- 
ing only the almost ideal part of the potential produc- 
ing the quadrupole trap potential and a second com- 
ponent, which contains mainly the effect of the trap 
imperfection under investigation. In our experiment we 
use a configuration with two diagonally opposite elec- 
trodes put on rf voltage J7if keeping the other pair at 
rf ground. The potential in this configuration is denoted 
by (0, C/rf)- It can be thought of as sum of the configura- 
tions (-C/rf/2,f7rf/2) and ([7^/2, f/rt/2). The configura- 
tion (— C/rf /2, C/rf /2) generates the quadrupole potential 
required for trapping and by symmetry has zero volt- 
age everywhere on the trap axis. Therefore, for the ideal 
symmetric trap, this configuration has no axial rf field 
component. The (J7rf/2, f7rf/2) configuration generates a 
fiat nonzero potential between the quadrupole electrodes 
with almost zero electric field in the neighborhood of the 
trap axis. It displays finite size effects and effects of the 
segmentation slits onto the on-axis electric field clearly. 

It can also be advantageous to apply this decompo- 
sition method when calculating the trap potential using 
the FEM. The precision achieved at a given mesh point 
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density increases with decreasing electric field. Thus, ef- 
fects which appear mainly in the ([/rf/2, U-ci/2) configu- 
ration, which does not contain the strong radial trapping 
field gradients, can be calculated with higher accuracy. 
The dependence of the accuracy of the solution on the 
electric field is accounted for by the definition of the 
residual error estimation function [49 , which plays an 
important role in approximating iterative linear solver 
methods employed in the solution of FEM problems as 
well as in adaptive mesh refinement routines. 

We estimate the accuracy of our calculations to re- 
solve the axial component of the rf field on the trap 
axis on a level well below 90 V/m at an rf amplitude 
C/rf of 1500 V. Regarding the large radial field gradients 
of 2.3 X 10® V/m^ obtained at this rf amplitude, this 
is demanding as the mesh elements in the region of the 
trap axis typically have a size in the range of 0.05 mm to 
0.1 mm and hence, fields on their surfaces are on average 
orders of magnitude larger. 

To reach the required accuracy we followed the mesh 
generation strategy with adaptive mesh refinement steps 



described in Section 4.1 until a convergence at the low 
percent level is reached. Where applicable, symmetry 
planes are used in the FEM model to reduce the size 
of the problem and hence required computing power. 
For the trap model of section 4.2 and 4.3, the linear de- 
composition into the configurations (— t/rf/2, Uii/2) and 
(C/rf/2, ?7rf/2) could be exploited and thus, higher ac- 
curacy reached at the same use of computing resources. 
For these calculations we obtained accuracies better than 
10 V/m for the axial component of the rf field on the trap 
axis. 



Dielectric components required in the actual imple- 
mentation of the trap were not taken into account in 
our field calculations, since they are either far away from 
the trap region or their surface is almost entirely coated 
with conducting electrode material. The electrode sur- 
faces are represented in the geometry model of the trap 
simulation by an assembly of rectangular boxes of 5 mm 
depth and of variable width 1^ and thickness as then 
stated in the text. An example of a geometric model 
used in a calculation for a trap made of three segments 
is schematically drawn in figure [3j 





4-.1 Technical details of the finite element calculations 

For the calculation of the rf trap potential, it is suffi- 
cient to consider the rf field at a fixed phase, such that 
the field calculation reduces to an electrostatic problem. 
Dynamic effects due to the electronic lead design that 
can produce different phase shifts for the rf field at dif- 
ferent points in the trap have been considered with a 
separate rf electronics software and will be discussed for 
the prototype trap [25|. 

The calculations of the rf field presented in this work 
consist in finding the solution to the electrostatic Dirich- 
let problem, in which the electrode surfaces of the trap 
are the boundaries kept at given electric potentials. We 
use a commercially available software (COMSOL Mul- 
tiphysics 3.5) for finite element analysis. In our FEM 
model we surround the trap electrode assembly with a 
wide rectangular boundary box, which is set to ground 
potential. The sides of the grounded box are typically 
at a distance of 1-2 cm from the trap electrode surfaces. 
We verified for different box sizes from 7 to 30 mm that 
this distance was large enough to have no influence on 
the obtained results. 



Figure 3: (a) Schematic drawing of a geometry model 
with three electrode segments used for FEM trap calcu- 
lations. The trap electrodes, including the grounded ends 
of the trap slit, are displayed as light blue rectangular 
boxes. Dedicated domains for mesh generation control 
are drawn in light to dark grey. The outmost boundary 
box set to ground potential is not shown, (b) Top view 
showing electrodes and domains of the top layer with 
grounded end electrodes. 



To generate the mesh we used the routines provided 
by the FEM software. The mesh elements are of tetrahe- 
dral shape. Quadratic Lagrange element functions were 
used in our calculations, with the exception of the cal- 
culations presented in Fig. [TT] as here the use of quartic 
Lagrange element functions led to a slight improvement 
in the quality of the solution. As control parameters for 
the mesh generation we mainly used the maximum mesh 
point separation and its growth rate in transitions from 
regions where a higher mesh point density is required 
to regions where it can be lower. In order to have an 
independent control over the mesh generation in differ- 
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ent regions of the geometry, we define various domains 
where dedicated mesh generation parameters can be ap- 
phed (see Fig. [s]). 

An initial mesh is generated with the maximum mesh 
point separation set equal to a fraction of the small- 
est distance in the surrounding local geometry and the 
element growth rate to values between 1.1 and 2. We 
changed mesh control parameters in the special mesh 
domains we had defined and compared the results of 
the calculation. We found that choosing small element 
growth rates in the range from 1.1 to 1.5, i.e. keeping 
mesh point density gradients small, was important to 
obtain convergence of the result with increasing mesh 
point density. Starting with a mesh, improved by the 
convergence analysis with manually changing mesh pa- 
rameters, we additionally apply the adaptive mesh re- 
finement routine of COMSOL. During the solution of the 
problem this routine detects the regions where improve- 
ments to the mesh are most effective in increasing the 
accuracy of the calculation. It accordingly improves the 
mesh quality and recalculates the solution. Typically af- 
ter two or three adaptive mesh refinement steps we found 
that, regarding the axial component of the field on the 
trap axis, the relative change between two consecutive 
adaptive mesh refinement steps had decreased to a level 
on the order of 10~^. 

The number of degrees of freedom or unknowns of the 
linear system of equations for our problem was typically 
in the range between 10^ and 8 x 10®. Only approximat- 
ing iterative solver methods are practicable for problems 
of this size and the computer memory at our disposal. In 
our case mostly an iterative conjugate gradient method 
with algebraic multigrid preconditioning and an iterative 
geometric multigrid method were used in combination 
with the direct solvers PARDISO and UMFPACK. As 
computer hardware we used two desktop computers, of 
which one has 4 gigabytes and the other has 16 gigabytes 
of main memory. 

4-8 Effect of the finite length of quadrupole electrodes 

For an electronically symmetric trap configuration 
{—Uy-f/2, Uy-f/2) the electrostatic potential along the trap 
axis is 0. Therefore, the electric field lines show no axial 
field component along the entire trap axis. For an highly 
segmented rf trap, where trapping electrodes are used si- 
multaneously for dc compensation, it is technically very 
challenging to deliver multiple independent rf voltages 
to the various electrode segments with identical ampli- 
tude and phase. Therefore, for the scalable trap design 
we chose the (0, U^f) configuration, where one rf trap 
electrode is kept on rf ground, facilitating independent 
adding of dc voltages. In this configuration the electro- 
static potential on the trap axis is roughly and 
electric field gradients along the axis are observed when 
terminating the rf trap electrodes to a finite size. In or- 
der to keep this axial component of the rf field amplitude 



below a given value, we adjust the minimum length of 
the trapping electrodes. For our linear Paul trap with a 
1 mm X 1 mm electrode spacing, figure |4]ja) shows the 
axial component of the rf field along the trap axis plot- 
ted over the axial coordinate z for various lengths of 
the quadrupole electrodes. Here, the electrode thickness 
= 0.2 mm and, for simplicity, an electrode geometry 
without notches was used. Our trap design also includes 
grounded electrodes at the ends of the quadrupole elec- 
trodes perpendicular to the trap axis behind a 1 mm 
insulating slit, as shown in figure |3] The purpose of the 
grounded end electrodes is to provide a defined electro- 
static potential between the ends of the quadrupole elec- 
trodes of the trap, where otherwise the insulator mate- 
rial of the wafer would be exposed. For the calculations 
shown in figure |4]ja) the grounded end electrodes are 
included in the trap geometry. For comparison, for an 
electrode length of 30 mm both cases are shown, with 
grounded and without grounded electrodes (free ends 
curve). 

Requiring an axial rf field component i?if,z of mag- 
nitude smaller than 90 V/m, we can find the width of 
the range on the trap axis usable for trapping of clock 
ions. In figure|4](b) the width of this usable range is plot- 
ted as a function of the length of the quadrupole elec- 
trodes. When we segment the quadrupole electrodes us- 
ing a segment length of 1 mm, for example, with a total 
length of the quadrupole electrodes of Lt — 30 mm and 
an electrode thickness of t^ = 0.2 mm, our design can 
accommodate a linear array of six usable ion traps for 
the clock application. The number of usable traps de- 
creases rapidly when the total electrode length is short- 
ened. With Lt — 2 cm and te = 0.2 mm the number 
of usable ion traps is already reduced to two. For an 
electrode thickness of 0.4 mm the corresponding usable 
range for clock application is given in figureQb) as well. 
In comparison to t^ = 0.2 mm, we obtain a moderately 
increased number of suitable traps: eight with an elec- 
trode length of 30 mm and three with a 20 mm electrode 
length. 

With a view to further miniaturization of trap assem- 
blies and further increase of the number of trap sites, it 
would be desirable to find a way to avoid long trapping 
electrodes. By introducing an additional rf ground elec- 
trode above and below the rf electrodes the effect of finite 
length of the trap electrodes can be decreased substan- 
tially, such that more trap segments can be used with the 
same or an even smaller length of trap electrodes. En- 
closing the linear quadrupole electrodes by extra ground 
electrode layers from above and below forces equipoten- 
tial surfaces of the rf potential to flatten out over a wide 
range of the trapping area, as shown in figure [5] (right) . 
With electric field lines being orthogonal on equipoten- 
tial surfaces, this results in a smaller gradient of -Eif.z 
over much of the z axis, with a steeper increase only 
close to the ends of the trap array. However, there are 
several drawbacks of the additional rf ground electrode 
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Figure 4: (a) Axial field compo- 
nent of rf field along the trap 
axis for different length Lt, cal- 
culated for an electrode thickness 
te = 0.2 mm. (b) Width of range 
along trap axis, where the axial 
component i?rf,z is of magnitude 
smaller than 90 V/m as a func- 
tion of the electrode length. Re- 
sults for an electrode thickness 
tf. = 0.2 mm and for — 0.4 mm 
are shown. Solid and dotted lines 
are to guide the eye. 
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Figure 5: False-color representation of rf equipotential 
surfaces of a linear quadrupole trap, (left) RF potential 
of a two-layer trap; (right) rf potential with additional 
rf ground electrode layers ontop of rf electrodes. 



layers, which are more or less important depending on 
the choice of electrode thickness and electrode layer spac- 
ing: the capacitance of the trap is increased; the numer- 
ical aperture for optical detection can become limited 
by the electrode assembly; increased risk of electrical 
breakdown; and finally, the axial rf field component on 
the trap axis resulting from notches in the electrodes is 
considerably intensified. As operating a 30 mm long ion 
trap in our setup brings no further limitation for the mo- 
ment, we did not include this option in our current trap 
design. Additional electrode layers are chosen only for 
extra compensation on top of the rf ground electrodes 
as shown in figure [2] that can then be put substantially 
closer. 

It is worthwhile mentioning that -Erf,z plotted in fig- 
ure |4]ja) can be well approximated in the central region 
on the z axis using the expression 



which gives the axial component of the field of a line 
charge of length L in the distance r from the z axis and 



with unit line charge density. For z values closer to the 
ends of the quadrupole electrodes the calculated field 
component -Erf.z deviates strongly from this expression 
as here the approximation using a constant line charge 
density can obviously not be made. Using a least square 
fit of this expression to the FEM calculation results for 
the axial electric field we obtain values for L, r and 
a proportionality factor C for different lengths of the 
quadrupole electrodes. In the calculations presented in 
the following paragraphs it was helpful to use a smaller 
length of the trap electrodes in order to keep the size 
of the problem smaller. To better isolate the effects un- 
der investigation we subtracted the finite length effect in 
form of the mentioned fit function from the FEM results. 

4- .3 Effect of insulation slits in dc electrodes 

In our design the rf ground electrodes also serve to apply 
dc voltages for axial confinement of the ions and to gen- 
erate electric compensation fields in axial and in one ra- 
dial direction. The electronic insulation of the electrodes 
requires segmentation slits. Here a notch is machined 
into the underlying insulator material and a conductive 
layer is applied around the teeth structure to cover in- 
sulation material in direct view of the ions, still leaving 
each tooth electrically insulated from its neighbors. 

Electrode surfaces of the trap assembly are equipo- 
tential surfaces of the electrostatic problem. The intro- 
duction of insulation slits into the electrode geometry 
causes a receding of equipotential surfaces in their vicin- 
ity. This receding is of course more prominent closer to 
the electrode. Further away, it is increasingly smoothed 
out and the curvature it causes in the equipotential sur- 
faces is reduced. Still, in the neighborhood of the trap 
axis a net effect of insulation slits in the electrodes is no- 
ticeable. Implementing corresponding symmetric notches 
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Figure 6: (a) Axial field com- 
ponent of rf field on the trap 
axis for different widths Wn of 
the notches in the quadrupole 
electrodes. The electrodes are 
notched at z = 0. The hori- 
zontal lines indicate the targeted 
maximum amplitude of 90 V/m. 
(b) Maximum magnitude of the 
axial component of the rf field 
|£'rf,z| as a function of w„. Lines 
between data points are to guide 
the eye. 
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in the rf electrodes can partially compensate for this 
effect. This can be understood when decomposing the 
corresponding electrostatic problem linearly into the two 
configurations (-C/rf/2, t/rf/2) and (C/rf/2, Urf/2). For the 
contribution of the (— C/i.f/2, [/i.f/2) configuration it is 
clear that keeping an exact geometric symmetry between 
the rf and rf ground electrodes leads to a zero potential 
all along the trap axis. On the trap axis, the effect of 
notches in the rf ground electrodes is then exactly can- 
celled out by the effect of the corresponding notches in 
the rf electrodes. However, even if the symmetry of the 
geometry is intact, there will be a residual magnitude in 
the axial component of the rf field which stems from the 
contribution of the (J7rf/2, C/rf/2) configuration and can 
not be compensated by symmetry. The potential on the 
trap axis in this configuration is a certain amount smaller 
than Urf/2, depending on how open the electrode struc- 
ture is. For trap electrodes, for example, being further 
apart from each other, the electrostatic potential on the 
trap axis is smaller. This seems to be exactly the effect of 
the notches in the electrodes, the equipotential surfaces 
are receded closer to the notches, leading to a minimum 
in the rf potential on the trap axis and resulting in a 
dispersion-shaped behavior of the axial component of 
the rf field. Figure |6] shows the calculation of this effect 
for different widths Wn of the notches in the electrodes. 
In this calculation a thickness tg of 0.2 mm was assumed 
for the electrode blades. To highlight the effect of a sin- 
gle segmentation, a trap geometry with only one notch 
per electrode blade was simulated. All four notches of 
the quadrupole electrodes are located at z = 0. 

We also calculated the effect for different electrode 
thickness. The maximum magnitude of the axial compo- 
nent of the rf field l-Brf.zl obtained for different electrode 
thickness is given in figure [7| for a fixed notch size 
Wn = 0.15 mm. Using thicker electrodes is obviously ad- 
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Figure 7: Maximum magnitude of the axial component of 
the rf field ji^rf z| as a function of the electrode thickness 
ig for notches of size w„ = 0.15 mm. Lines between data 
points are to guide the eye. 



vantageous regarding the contribution of the notches to 
the axial component of the rf field. This can be an in- 
teresting option when the electrodes are made using a 
technique different from the one we chose. Both laser 
cutting of thin wafers and the subsequent metallization 
are more difficult for thicker structures and higher aspect 
ratios. 

For a segment length of a few mm and smaller, the 
effects of the notches from both ends of a segment can- 
cel partially. The behavior of the resulting axial rf field 
in the middle of the segment varies with the segment 
length, since it is composed of contributions from both 
ends of the segment, which have opposite sign. Figurejsjja) 
shows the simulation of a trap with three segments, thus 
two discontinuities in the electrode blades. The inner seg- 
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Figure 8: (a) Axial component of the rf field i?rf,z on 
axis as a function of the axial coordinate z in a single 
trap for a segment length Is = 1 mm and for a segment 
length Zg = 2 mm with the width of the notches being 
0.15 mm. (b) Axial component of the rf field i?rf.z on 
axis as a function of the axial coordinate in an array 
of traps with 5 inner segments of length Ig ^ I mm 
separated by 0.15 mm wide notches. The axial location 
of the segments is indicated by dotted lines. 



ment length is Zg = 1 mm and = 2 mm, respectively, 
with the width of the notches kept at 0.15 mm and elec- 
trode thickness at 0.2 mm. Here, the center of the inner 
segment is at z = 0. The curve for the 1 mm long seg- 
ment shows, that the effects of the notches from both 
sides of the segment are partially cancelling each other. 
This can also be seen in the calculation for a trap with 
five inner segments of 1 mm length shown in figure [sj^b). 
We can conclude, that if the symmetry of the electrode 
construction is given, the remaining axial rf field compo- 
nent produced by the notches in the electrodes is small 
enough to allow trapping of ion chains of five to ten 
ions for clock applications. Even if a larger width of the 
notches is chosen, close to the center of each trap seg- 
ment the axial rf field component has a zero crossing 
resulting in a range with l-Brf.zl < 90 V/m large enough 
to accommodate a chain of ions. 



4.4 Angular misalignment of quadrupole electrodes 

Due to the precision achieved in the process of laser ma- 
chining of 1 or 2 /im, we can assume that electrodes lo- 
cated on one wafer can be aligned nearly perfectly with 
respect to each other. Thus, we only need to consider 
angular misalignment of the two quadrupole electrode 
pairs occurring during the trap assembly. 

First, we consider a possible rotation angle a between 
the two wafers of the quadrupole trap electrodes around 
the axis normal to the wafer plane, corresponding to the 
X axis in figure [2] Figure [9f^a) shows the axial compo- 
nent of the rf field calculated for an electrode thickness 
of 0.2 mm for different misalignment angles a using a 
(0, C/rf ) configuration. Here, the trap axis 2' is rotated by 
the angle a/2. Our calculations show that this angular 
misalignment produces an rf field component along the 
trap axis, which is proportional to the angle of rotation of 
the two electrode layers with respect to each other. This 
can be seen from figure ^h) which shows a plot of -Erf ,z' 
as a function of the rotation angle. When repeating the 
calculation for an electrode blade thickness of 0.4 mm 
the results for the axial rf field component show no sig- 
nificant difference compared to the results obtained with 
0.2 mm electrode thickness. We deduce from figure |9|^b) 
an axial rf field of 6.2 x 10^ V/m per mrad. In order 
to keep this axial rf field component well below 90 V/m 
we have to ensure that the angular misalignment of the 
two wafers for rotation around the x-axis is less than 
0.14 mrad. 

Second, we consider an angular misalignment of the 
two electrode wafers with respect to each other around 
the axis of rotation corresponding to the y axis (see 
fig. [2|. An angular misalignment around this axis can 
occur during the assembly of the electrode wafers, when, 
for example, the spacer height is not identical, or glue is 
creeping between spacer and electrode blade. Such a mis- 
alignment leads in first order to a linearly changing dis- 
tance between the electrode wafers along the trap axis. 
This obviously leads to a corresponding change in the 
radial rf field gradients and consequently of the radial 
trap frequencies for the secular motion and a small vari- 
ation of the on-axis potential is to be expected as well. In 
our calculations we used misalignment angles in a range 
from 1 mrad to 10 mrad. However, within the estimated 
accuracy of the calculation, which is well below 90 V/m 
on the trap axis, we could not resolve any deviation from 
zero in the axial rf field component due to angular mis- 
alignment of this type. 



4-5 Tolerances on width and position of notches in 
electrodes 

In the process of machining notches into the electrode 
blades inaccuracies can occur in both the positioning of 
the notches as well as in their width. The subsequent 
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Figure 9: (a) Axial component 
of the rf field i?rf,z' on the trap 
axis {z') as a function of the axial 
coordinate for different misalign- 
ment angles a. Here the trap axis 
(z') is rotated by the angles a/2 
and —a/2 to the axes defined 
by the upper and lower electrode 
layer, respectively, (b) Magni- 
tude of the axial component of 
the rf field i?rf,z' on the trap axis 
as a function of the misalignment 
angle. Data points are fitted with 
a linear fit to extrapolate i?rf ,z' to 
small angles. 
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process of coating the insulator with the electrode ma- 
terial can also contribute to further inaccuracies. 

First, we consider the effect of a small change in the 
width of one of the four notches separating two segments 
of the linear Paul trap. We calculated the resulting rf 
field when one notch is made wider in one rf carrying 
electrode. Our calculations show what can be expected 
from such a change in the geometry: an additional reced- 
ing of the equipotential surfaces away from the widened 
notch occurs, centered at the axial coordinate of the 
notch. As a result the axial rf field component displays a 
dispersion-shaped behavior along the z-axis at the edge 
of the trap segment in excess of the magnitude found for 
a symmetric electrode geometry as shown in figure [6f^a) . 
The radial components of the rf field as well display on 
the z-axis an increase from zero which peaks at the z 
location of the widened notch. The feature in the radial 
field components is of similar width and amplitude as 
the one in the axial component. Due to the strong gra- 
dient in the radial rf field components the excess field 
produced by the widened notch leads to a small bending 
of the actual trap axis away from the z-axis by a distance 
on the order of 100 nm. 

Figure [lO{^a) shows the axial rf field component as 
a function of the axial coordinate obtained for differ- 
ent positive deviations Aw in the width of one notch 
from the width of the notches in the other three elec- 
trodes, which are kept at Wn = 0.1 mm. The depen- 
dence of the amplitude of the dispersion-shaped feature 
in the axial rf field component on the deviation in the 



width of one notch is shown in figure 10 b). The line 



in the graph resulting from a least square fit suggests, 
that the dependence can be regarded to be linear in this 
range. Assuming a segment length of Is — 1 mm, we de- 
duce a maximum tolerance for the notch width in the 
rf carrying electrode of Aw = 3 /xm in order to have 



|£'rf,z| £ 90 V/m in the middle of the segment. For a 
segment length of 2 mm the middle of the trap segment 
is further away from the segment border and the region 
where -Erf,z peaks, such that notches can be cut to less 
stringent tolerances of Aw = 11 /zm. 

By widening a notch in an rf ground electrode instead 
of in an rf carrying electrode we find in our calculations 
a very similar effect in the axial component of the rf field 
on the trap axis, albeit smaller in amplitude, by roughly 
20% for the here considered geometry, and of opposite 
sign. 

We also varied the design width Wn of the electrode 
notches in the range from 0.1 mm to 0.2 mm and ob- 
served a linear increase of the amplitude of the effect 
of a widened notch with increasing w„. The amplitude 
of the effect was increased 1.8 times for a design notch 
width of 0.2 mm as compared to a design notch width 
of 0.1 mm. 

The calculations in figure [TO] have been carried out 
for an electrode thickness of ~ 0.4 mm. When we 
repeated the calculations for te = 0.2 mm we found only 
a small increase in the magnitude of i?rf,z by a factor of 
1.07 for a design notch width of Wn — 0.1 mm and by a 
factor of 1.15 for w^ = 0.2 mm. 

Second, we consider an inaccuracy Az in the position 
of one notch in an otherwise perfectly symmetric elec- 
trode assembly. Figure 11 shows the effect on -Erf.z when 
one notch is displaced in positive z direction. On the 
axial component of the rf field the effect is asymmet- 
ric with a maximum located around the center of the 
notches followed by a minimum of smaller magnitude 
located around z = 0.7 mm. The magnitude of the ad- 
ditional axial rf field produced by a position inaccuracy 
is considerably smaller than that found for inaccuracies 
in the width of a notch. Since the maximum of l-Erf z| 
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Figure 10: (a) Axial component 
of the rf field E'lf.z on the trap 
axis as a function of the axial co- 
ordinate for different deviations 
Aw in the width of one notch 
in an rf carrying electrode. All 
other notches are kept at Wn — 
0.1 mm. The notches are cen- 
tered at z = 0. (b) Maximum 
amplitude of i?rf.z along the trap 
axis as a function of Aw. Data 
points are fitted to show their lin- 
ear dependence in this range. 
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Figure 11: (a) Axial component 
of the rf field E'rf.z on the trap 
axis as a function of the axial co- 
ordinate for different deviations 
Az in the position of one notch. 
In the geometry used for the cal- 
culations Wn = 0.1 mm, except 
for the curve mentioned in the 
legend, where w^ = 0.2 mm. 
(b) Maximum of |i?rf,z| as a func- 
tion of Az. 
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is also localized on the border of the segment the trap 
region in the center of the segment is even less affected. 



For the calculations shown in figure 11 a.) the notch 
width Wn =0.1 mm has been used with exception of 
one curve, for which the chosen notch width of 0.2 mm 
is pointed out in the legend of the graph. We observed a 
rapid increase of the maximum of l-Erf.zl as a function of 
the notch width, the displacement of one notch with re- 
spect to the others being kept constant. We want to point 
out that the residual axial rf field amplitude occurring 
in the perfectly symmetric case also grows rapidly with 
increasing width of the notches (see fig. |6| and can con- 
tribute significantly here. The curves in figure 11 [a) have 



been calculated using an electrode thickness of 0.4 mm. 
When we used an electrode thickness of 0.2 mm instead, 
we did not find a significant increase in the effect. 



Targeting a positioning accuracy Az of better than 
10 /im when aligning the different electrode layers on 
top of each other, the axial rf field contribution due to 
position inaccuracies will not become critical for a notch 
width of up to 0.2 mm. 



4.6 Translational misalignment of quadrupole electrodes 

We looked into the misalignment of the trap assembly 
in which one of the electrode wafers is translated with 
respect to the other parallel to the trap axis. The trans- 
lation destroys the symmetry of the geometry at the end 
of the trap segments, where the notches in the electrodes 
of one wafer become displaced with respect to the corre- 
sponding notches in the other wafer. Considering such a 
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pair of notches of which each is in another wafer, we ex- 
pect that the displacement resuhs in an extra axial com- 
ponent of the rf field halfway between the two notches. 
In our FEM calculations the extra axial rf field contri- 
bution was evident between the wafers in the vicinity of 
the notches. 

The calculations also show that the extra axial field 
decreases in magnitude towards the trap axis to finally 
vanish right on the axis. Since for the pair of notches 
situated on the opposite side of the trap axis the con- 
figuration of rf and rf ground electrodes is reversed on 
the two wafers, the additional axial rf field component 
has opposite sign on both sides of the trap axis. As a re- 
sult the effects of the translation from two opposite pairs 
of notches cancel each other on the trap axis. From our 
calculation using w„ = 0.1 mm and = 0.4 mm we con- 
clude that even at a distance of 0.05 mm from the trap 
axis and at z coordinates corresponding to the location 
of electrode notches, an axial translation of 0.05 mm of 
the wafers with respect to each other results in a con- 
tribution to the axial component of the rf field smaller 
than 90 V/m. In the middle of the segment where the 
ions are trapped the effect is even smaller. 



4.7 Introduction of extra electrode layers for dc field 
compensation 

In order to compensate for dc stray fields in radial di- 
rection, which would displace trapped ions from the null 
field line of the rf field, the trap design has to offer a way 
to apply independent dc compensation fields in the ra- 
dial plane in each trap segment. Two degrees of freedom 
are required to compensate for any electric field in the 
radial plane. Using a small voltage difference between 
the rf ground electrodes of the trap segment generates a 
dc electric field oriented in the (x,y) plane nearly along 
the (1,1) direction. In order to generate a second dc elec- 
tric field with a strong component perpendicular to this 
axis, we considered two alternative designs. 

In a first design, schematically drawn in figure [12] ^ a), 
we added segmented compensation electrodes behind a 
1 mm wide rf electrode with an insulation gap of 0.2 mm 
between rf and compensation electrodes. This design is 
compact but bears some strong shielding of the com- 
pensation electrodes by the trap electrodes. The FEM 
calculation of the dc field generated in this geometry 
with an electrode thickness of 0.4 mm, a segment length 
of 1 mm and a notch width of 0.1 mm gives a field of 
8 V/m in the trap region of the segment using -1-1 V 
and -1 V on the extra compensation electrodes of the 
segment. Also in the next and second to next neighbor- 
ing segments considerably large dc fields of 4 V/m and 
1.4 V/m are found. 

A second design consists in expanding the two-layer 
quadrupole trap assembly by two additional layers with 



a) 



□ rf ground / dc 2 

□ rf ground / dc 1 



b) 



I rf ground / dc 1 
I rf ground / dc 2 
I rf 



Figure 12: Alternative electrode designs to obtain two 
degrees of freedom for dc field compensation in the radial 
(x, y) plane. 



gives a schematic view of this design. The additional lay- 
ers are similar to the other two trap layers only that the 
rf electrode side is missing. Each extra layer is placed 
with a spacing of 0.25 mm on top and precisely aligned 
with a segmented rf ground electrode. The geometry of 
the extra compensation layers is kept identical to the in- 
ner trap electrodes regarding notch width and electrode 
thickness. The FEM calculation of the dc field generated 
by this geometry using an electrode thickness of 0.4 mm, 
a segment length of 1 mm, and a notch width of 0.1 mm 
gives a stronger field value in the trap region of 60 V/m 
for the voltages of -1-1 V and -1 V on the two extra com- 
pensation electrodes. The direction of the generated field 
vector points mainly along the a;-axis with a component 
of 3.7 V/m along the y-axis. The corresponding field val- 
ues obtained in the next and second to next neighboring 
segments are 19 V/m and 2.5 V/m, respectively. 

The second design has a clear advantage over the first 
design regarding the larger field. The field is also well 
localized on the trap segment with only a 4% contribu- 
tion to the dc field in the neighboring traps. In addition, 
considering the relatively large rf voltage amplitude we 
envisage to apply, the second design is also less prone to 
rf breakdown and surface fiashover than the first. The 
slightly larger layer stack reduces the solid angle for op- 
tical detection to an acceptable level. 

Since we opt for the second design we investigate the 
effect of the extra compensation electrode layers on the 
rf field. It is important to examine how the additional 
compensation electrodes affect the axial component of 
the rf field produced by the electrode notches. Consid- 
ering the effect on a single segment border we calculate 
the rf field for an electrode geometry with only two long 
segments, thus one central discontinuity in the electrode 
blades. Adding the extra compensation electrode lay- 
ers to the trap geometry clearly increases the amplitude 
of the axial rf field component caused by the notches. 
Compared to calculations for the same geometry with- 
out extra compensation electrodes as described in sec- 
(figure [6]) , we find that the amplitude of i^rf 



4.3 



dedicated segmented compensation electrodes. Fig. 12 ^b) 



tion 

is increased by a factor of 3.6 for an electrode thick- 
ness of 0.2 mm and a notch width of 0.2 mm, and by 
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a factor of 4.6 for an electrode thickness of 0.4 mm and 
notch widths of 0.1 mm as well as 0.15 mm. Although 
the extra compensation layers lead to a considerable in- 
crease in the axial rf field component on the segment 
borders, the axial rf field component can still be small 
enough in the region in the middle of a trap segment, 
because here again the effects stemming from neighbor- 
ing segment borders partially cancel one another out. An 
example for this is given in figure [13] for three different 
design parameters: electrode thickness, notch width and 
segment length. Here, a three segment trap was simu- 
lated with the middle segment of length Is surrounded 
by two long outer segments. In the calculations for all 
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Figure 13: (a) Axial component of the rf field i?rf,z on 
the trap axis as a function of the axial coordinate z for 
a 3-segment trap design with extra compensation elec- 
trode layers. The center of the middle trap segment is 
located at 2; = 0. The curves are calculations for differ- 
ent geometrical design parameters: electrode thickness 
^e, notch width and segment length 1^ given in mm 
in the legend. 



three parameter sets shown in figure [13] an axial region 
of at least 0.12 mm length with l-Erf.zl < 90 V/m can be 
found, which is suitable for clock spectroscopy. 



5 Conclusion and outlook 

We presented an electrode design of a segmented linear 
Paul trap optimized for optical clock application and fre- 
quency standards based on the spectroscopy of Coulomb 
crystals. Using segmented electrodes the design supports 
the implementation of linear arrays of traps, where in 
each trap segment electrodes are provided for the com- 
pensation of dc electric fields in all three directions of 
space. The trap assembly still offers a good optical ac- 
cess for laser beams and fluorescence detection. 



In this design geometry we propose to trap arrays 
of ^^^In+ ions sympathetically cooled by ^''^Yb"'". We 
showed that with such an optical frequency standard 
of improved short-term stability a fractional long-term 
stability of 10~^^ can be reached with systematic fre- 
quency shifts at the mHz level. Integration times for the 
frequency measurement can be shortened by more than 
an order of magnitude and requirements for clock laser 
stability relaxed significantly. 

Of great importance for high precision spectroscopy 
of a larger sample of ions is the minimization of excess 
micromotion of the ions, which can contribute consid- 
erably to the second-order Doppler shift of the clock 
transition. In our design we minimized the critical ax- 
ial component of the rf field along the trap axis and 
studied tolerances for the manufacturing process of the 
trap. We identified and quantified the three most critical 
causes for excessive magnitudes of the axial rf field com- 
ponent: the finite length of the quadrupole electrodes; 
the tolerance on the width of the notches separating elec- 
trode segments; and the tolerance on the alignment angle 
of opposing electrode wafers (torsion around the x axis 
in our coordinate system). For our choice of geometry 
with an 1 mm x 1 mm spacing between the quadrupole 
electrodes, and requiring the axial component of the rf 
field on the trap axis to not limit the targeted fractional 
uncertainty of the frequency standard of 10^^^, we can 
summarize our results on these critical effects as follows. 
Using a trap segment length of 1 mm and electrode thick- 
ness of 0.2 mm a linear array of six usable traps with 
at least 10 ions each is obtained with a total electrode 
length of 3 cm, whereas with a total electrode length of 
2 cm our design only supports two usable traps. Choos- 
ing the width of the notches that are separating trap seg- 
ments to be 0.1 mm and the segment length to be 1 mm, 
the tolerance on the width of a notch in the electrodes 
is 3 fiia. The tolerance on the angular misalignment of 
the two quadrupole electrode layers around the direction 
orthogonal to the wafers is less than 0.14 mrad. 

The effect of other tolerances in the trap construc- 
tion, like the position deviation of a single notch in the 
electrodes, the translational displacement along the trap 
axis of one electrode layer with respect to the other, 
and the angular misalignment of the electrode layers 
around the y axis, have also been investigated and quan- 
tified. We showed that their effect on the axial compo- 
nent of the rf field on the trap axis is either negligibly 
small or can be easily avoided by applying state-of-the- 
art tolerances in the manufacturing process. Further, we 
point out that depending on the choice of segment length 
boundary effects of the segmentations can cancel each 
other out, reducing the on-axis micromotion. 

Based on our FEM calculations a prototype trap 
with two trapping segments made of Rogers4350B"'"'^ 
has been constructed and is currently under test. First 
measurements of micromotion with ^^^Yb"*" ions show a 
promising low value in agreement with our simulations. 
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For a scalable high precision trap we are developing a 
laser cutting and metalhzation process for AIN wafers, 
in which we plan to investigate an array of 6 x 10 ions 
for an optical frequency standard. 
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